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Universality class of fiber bundles with strong heterogeneities 
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Abstract. - We study the effect of strong heterogeneities on the fracture of disordered materials 
using a fiber bundle model. The bundle is composed of two subsets of fibers, i.e. a fraction 
< a < 1 of fibers is unbreakable, while the remaining 1 — a fraction is characterized by a 
distribution of breaking thresholds. Assuming global load sharing, we show analytically that 
there exists a critical fraction of the components ct c which separates two qualitatively different 
regimes of the system: below a c the burst size distribution is a power law with the usual exponent 
t — 5/2, while above a c the exponent switches to a lower value r = 9/4 and a cutoff function 
occurs with a diverging characteristic size. Analyzing the macroscopic response of the system 
we demonstrate that the transition is conditioned to disorder distributions where the constitutive 
curve has a single maximum and an inflexion point defining a novel universality class of breakdown 
phenomena. 
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J Introduction. — Damage and fracture of materials 
occurring under various types of external loads is a very 
important scientific problem with an enormous technolog- 
ical impact. During the last two decades the application of 
statistical physics has revealed that heterogeneities of ma- 
terials' microstructure play a crucial role in fracture pro- 
cesses [1]. To capture the effect of disorder, recently sev- 
eral stochastic fracture models have been proposed such 
as the fiber bundle model (FBM) and lattice models of 
fuses or springs [1-7]. Based on these models, analytic cal- 
culations and computer simulations revealed that macro- 
scopic fracture of disordered materials shows interesting 
analogies with phase transitions and critical phenomena 
having several universal features independent of specific 
material details [1,4-6,8-10]. It has been found that un- 
der a slowly increasing external load macroscopic failure 
is preceded by a bursting activity due to the cascading 
nature of local breakings [3,4]. Since the bursts can be 
recorded experimentally by the acoustic emission tech- 
nique, these precursors addressed the possibility of fore- 
casting the imminent failure event [11-14]. The size dis- 
tribution of bursts was proven to be a power law with an 
exponent which is universal for a broad class of disorder 
distributions [3,4]. Recently, the robustness of the uni- 
versality class has been tested by mixing different types 
of disorder distributions [17], and by introducing a gap 
into the domain of strength values [18]. However, relevant 
change of the burst size distribution was only obtained 
when introducing a finite lower threshold for the strength 



disorder. Increasing the threshold strength a crossover oc- 
curs from a power law of exponent 5/2 to another one with 
a lower exponent 3/2 [14-16]. Divakaran and Dutta have 
studied the critical behaviour of a Random Fiber Bun- 
dle Model with mixed uniform distribution of threshold 
strengths [18]. They have considered two uniform distri- 
butions separated by a gap. The approach developed in 
this Letter might be interpreted as the infinite gap limit 
of Divakaran's model. 

In the present paper we study the effect of strong het- 
erogeneities on the process of fracture based on a fiber 
bundle model. We assume that the system has two com- 
ponents one of which is characterized by a strength distri- 
bution, while the other one is unbreakable. Varying the 
fraction of the two components a under global load shar- 
ing conditions, we show analytically that the presence of 
unbreakable elements has a substantial effect on the frac- 
ture process of the system both on the micro- and macro- 
scales. Very interestingly, we find a critical fraction a c 
where a transition occurs between two qualitatively dif- 
ferent regimes: below the critical point a < a c the macro- 
scopic constitutive curve has a single maximum and the 
burst size distribution is a power law with the usual mean 
field exponent r = 5/2. However, above a c the macro- 
scopic response becomes monotonous and the burst expo- 
nent switches to a lower value r = 9/4 with a cutoff func- 
tion. Based on the analysis of the macroscopic response 
of the system, we show that the transition is conditioned 
to disorder distributions where the constitutive curve has 



R. C. Hidalgo et al. 



a single maximum and an inflection point defining a novel 
universality class of breakdown phenomena. 

Model. — We consider a set of N fibers which are 
loaded in parallel. Under an increasing external load u a 
the fibers have a linearly elastic response with a Young 
modulus E = 1 fixed for all the fibers. In order to cap- 
ture the large variation of disordered material properties, 
we assume that the bundle is composed of two subsets of 
fibers with strongly different breaking characteristics: A 
fraction < a < 1 of fibers is strong in the sense that they 
have an infinite load bearing capacity so that they never 
break. However, fibers of the remaining f — a fraction are 
weak and break when the load on them a exceeds a thresh- 
old value a\ hl i = 1, . . . , N w , where N w = (1 — a)N is the 
number of weak fibers. The strength disorder of weak 
fibers is characterized by the probability density p(ath) 
and distribution function P(o~th) = J ±h p(x)dx of the fail- 
ure thresholds. After a weak fiber breaks in the bundle, 
its load has to be overtaken by the remaining intact ones. 
For simplicity, we assume global load sharing (GLS) (also 
called equal load sharing) which means that all the intact 
fibers share the same load a, hence, no stress concentra- 
tion occurs around failed regions. Under these conditions 
the constitutive equation of the model can be written as 

a = (l-a)[l-P(a)}a + ao-, (1) 

where a is the external load acting on the sample and 
a denotes the load of single fibers which is related to the 
strain e of the system as a = Ee. The first term of Eq. 
(1) accounts for the load bearing capacity of the surviv- 
ing fraction of weak elements, and the second one repre- 
sents the stress carried by the unbreakable subset of the 
system. In the following calculations it is instructive to 
consider two different strength distributions for the weak 
fibers, namely, a uniform distribution between and 1 
and a Wcibull distribution will be used with the distribu- 
tion functions -P(er) = a and P(cr) = 1 — cxp [— (er/A) m ], 
respectively. 

The constitutive behavior of the system is presented in 
Fig. 1, for the two different disorder distributions. We 
recover the usual FBM solutions [2] in the limiting case 
of a = 0, when the bundle is only composed of weak 
fibers. Those solutions usually present a parabolic maxi- 
mum, which defines the critical deformation a c = Ee c and 
critical strength u^(u c ) of the system. For finite values of 
a, all the weak fibers break for large enough a so that 
the first term of Eq. (1) goes to zero while the unbreak- 
able fibers overtake the entire external load. Consequently, 
the constitutive curves in Fig. 1 tend asymptotically to a 
straight line with slope aE. It can be seen in Fig. 1 that for 
low values of a the local maximum of 0~ o {a) prevails but its 
position o~ c {a) and value tr^a) are monotonically increas- 
ing with a. It is interesting to note that there exists a well 
defined critical value of the fraction of the components a c 
above which a > a c the local maximum disappears and 
the constitutive curve becomes a monotonically increasing 
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Fig. 1: Constitutive behavior of the system for several values of 
a using a Weibull distribution with m — 2 and A = 1 (a), and 
uniformly distributed threshold values (b). The two regimes 
a < etc and a > a c are indicated by the grey areas, and the 
dashed lines with slope a show the asymptotic linear behavior 
of a (a). In the figures, the vertical straight lines represent 
the position of the inflexion point and the value of the largest 
breaking threshold for the Weibull (a) and uniform distribu- 
tions (6), respectively. 



function da /da > 0. The position of the maximum a c is 
obtained from the condition of extreme ^r 2 - 1 = 0, which 

d(T \<T C 5 

can be cast into the form 

-^— = P{a c ) + a cP {a c ). (2) 
I — a 

The above equation should be solved for a c as a function 
of a, then a c Q can be determined by substituting a c (a) 
into Eq. (1). Since the derivative of the constitutive curve 
o~ (o~) has a minimum in the inflexion point ai„, the right 
hand side of Eq. (2) has a maximum at a c = <ii n . It 
follows that Eq. (2) can only be solved for a c (a) until 
a < a c , where the critical fraction of strong fibers a c is 
defined as the solution of cr c (a c ) = <Ji n . It is important to 
emphasize that the location of the inflexion point erj„ does 
not depend on the value of a since the second derivative 
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of the constitutive equation a (a) reads as 



analytically in the form of an integral 



d 2 a n 



d 2 c 



-(1-a) [2p(a) + ap'(a)} 



(3) 



For the case of Weibull distributions the general solution 
a c (a) cannot be obtained in a closed form. However, one 
can still calculate analytically a c and determine a c for 
two parameter values a = and a = a c . The calculations 
result in a c (a = 0) = A(l/m) 1 / m , and a c (a c ) = a in , where 
the inflexion point is ai n = A [(1 + m) / m] 1 ^ r ™. The critical 
point a c is then obtained by substituting tTi n into Eq. (2) 
which yields a c = me -(i+™)/™/ [i + me -(i+m)/ m ] _ 

Note that the above arguments do not apply to the 
uniform distribution, since the constitutive curve a = 
[1 — (1 — a) a] a does not have an inflexion point (see Fig. 
1(6)). The position of the maximum of a Q (a) can be ob- 
tained analytically as a c (a) = 1/2(1 — a), which holds for 
a < a c with the critical value of the control parameter 
a c = 1/2. At a c the value of a c coincides with the up- 
per bound of strength values o™^ x = 1. The parabolic 
shape of the constitutive curve prevails even for a > a c 
but a (a) becomes linear at a = a™ x before reaching 
the maximum, so that the rest of the parabola cannot be 
realized. 

The presence of the critical point and the qualitatively 
different forms of a (a) below and above a c have a sub- 
stantial effect on the microscopic breaking of the system. 
Under stress controlled loading conditions the decreasing 
part of a (a) can not be accessed for a < a c . Contrary, 
a horizontal jump occurs giving rise to a large number 
of breakings in one step. For uniformly distributed fail- 
ure thresholds, this unstable avalanche is the last one that 
includes all the remaining weak fibers. In the Weibull 
case, however, the threshold values are distributed over 
an infinite domain, so that the jump is still followed by 
breaking events which disappear only asymptotically. For 
the detailed characterization of the microscopic breaking 
process, we analyze the size distribution of bursts of fiber 
breakings. 

Precursory activity. — Under stress controlled load- 
ing conditions, each fiber breaking is followed by the re- 
distribution of load over the intact elements. Assuming 
global load sharing the load is everywhere the same a in 
the system. When the external load is increased quasi- 
statically, i.e. a is increased to break only a single fiber, 
the subsequent load redistribution triggers an entire burst 
of breakings. In the simple FBM these local failure events 
result in fluctuating burst sizes A, with an increasing aver- 
age, as macroscopic failure is approached. The size distri- 
bution of the bursts is one of the most important charac- 
teristics of the microscopic fracture process which can be 
monitored experimentally by the acoustic emission tech- 
niques. It has been demonstrated that in FBM under GLS 
conditions [4] , the burst size distribution can be obtained 
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p(a)(l-a ff )a^- 1 e-^ A da, (4) 



where a a = ap(a)/[l — P(a)] is the average number of 
fibers which break as a consequence of a single fiber failure 
at the load a. It was shown in Refs. [4, 14, 15] that the 
distribution D(A) simplifies to a power law D(A) ~ A~~ r 
with the exponent r = 5/2 for a broad class of disorder 
distributions where the constitutive curve of the system 
has a single quadratic maximum. 

In the following we show analytically that in the pres- 
ence of unbreakable fibers, the avalanche statistics changes 
and a novel universality class of FBMs emerges. Slowly in- 
creasing the external load to break a single fiber, its failure 
stress a is equally redistributed over the intact fibers giv- 
ing rise to the load increment 



Sa 



N[l- P(a)] (1 - a)+aN' 



(5) 



It can be seen that the strong fibers reduce the load incre- 
ment Sa on the weak ones, since the load beared by the 
strong fibers does not contribute to breaking. The aver- 
age number of fibers a a which fail as a consequence of this 
increment Sa can be cast into the form 



a a = (1 — a)Np(a)Sa 



(1 — a)ap(a) 



a+(l-a) [1-P(a)}' 



(6) 



The size distribution of the resulting bursts can be ob- 
tained by substituting Eq. (6) into the general expression 
Eq. (4), where we have to analyze the behavior of the in- 
tegral 
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^ -A[o„-lnn„ 



da. 



(7) 



for different values of a. The upper integral limit corre- 
sponds to the location of the maximum in the constitu- 
tive curve a rn . For large A this integral is controlled by 
the maximum of the exponent. The extreme condition of 
ip = a a — In a a result in ip' = a^(l — — ) = 0, correspond- 
ing to a maximum at a a = 1. Below the critical point 
a < a c , for a a < 1, we can make the expansions 



and 



(8) 
(9) 



<i/j ~ 1 

Inserting Eqs. (8,9) into the expression of /(A) we get 

\a-a m )e-^^'^ 2 da. (10) 



1(A) ~ p(a m )a (J e 



Substituting 1(A) into the general equation Eq. (4) and 
taking the large A limit of the prefactor, the asymptotic 



R. C. Hidalgo et al. 



behavior of the burst size distribution can be cast in the 
form 

D(A) _ p{a m ) A _ 5/2 



N 



2na' 



-A- 



(11) 



which coincides with the known result of Refs. [4, 14, 15] in 
the limit a = 0. This derivation implies that the presence 
of a finite amount of unbreakable fibers does not change 
qualitatively the behavior of the system while the single 
quadratic maximum of the constitutive curve prevails a < 
a c . 

The situation drastically changes when we reach a c , 
since at this point the position of the maximum cr c (a c ) 
and of the inflexion point <Ji n of the constitutive curve co- 



incide with each other so that ^ 







= 0and *y 

hold. Above a c no maximum of the constitutive curve ex- 
ists da jda > 0. It can easily be shown that at a — a c 
the average number of failing fibers a a as a consequence 
of a single fiber breaking has the properties a (Tc = 1 and 
a' a | CTo = 0. In order to determine the asymptotic behavior 
of /(A) it is then necessary to carry out the Taylor ex- 
pansions of Eqs. (8,9) to the next order. In this case we 
obtain for a a and ip 



a a ~l + ^{a-a c )\ ,p ^ i + ^ {<7 - a c )\ (12) 

Inserting these expressions into Eq. (7), the integral can 
be cast into the form 



Following the usual procedure, we arrive at 

1" (I) „-A^-3/4 



1(A) 



2 5/2 3 7/4 a /,l/2 



(13) 
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which implies that at the critical point a c the avalanche 
size distribution changes to 
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Our derivation demonstrates that increasing a the behav- 
ior of the system changes both on the macro- and the 
micro-scales. We showed that while the quadratic maxi- 
mum of (J (&) prevails, i.e. below the critical point a c , the 
asymptotic behavior of the burst size distribution D(A) 
is controlled by the vicinity of the maximum resulting in 
a power law functional form -D(A) ~ A _T with an uni- 
versal exponent r = 5/2. However, at a c the constitutive 
curve becomes monotonically increasing da Q /da > and 
the avalanche statistics is dominated by the inflexion point 
of <x ((x), giving rise to a different value of the exponent 
t = 9/4. Varying the control parameter a, the exponent 
t suddenly switches between the two values 5/2 and 9/4 
when passing the critical point a c . Note that in the deriva- 
tion the only assumption we made is that the constitutive 
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Fig. 2: Non-normalized avalanche size distributions for uni- 
form (a, c) and Weibull distributions (b, d) varying a below 
and above a c . The straight lines in (a) and (6) represent the 
power laws obtained analytically Eqs. (11,15). Rescaling the 
two axis according to the scaling formula Eq. (16), a very good 
quality data collapse is obtained in (c) and (d). 



curve of the system has a single maximum and an inflexion 
point. It follows that the change of the exponent t of the 
avalanche size distribution can be observed for a large va- 
riety of disorder distributions defining a novel universality 
class of breakdown phenomena. This universality class is 
narrower than the one in which the power law behavior of 
D(A) emerges with the exponent r = 5/2. For instance, 
the Weibull distributions do present the above switching 
of exponents, however, the uniform distribution does not. 

We have carried out Monte Carlo simulations to validate 
the previous theoretical predictions. We have explored 
the quasi-static fracture process of our fiber bundle model 
using computer simulations of a system composed of N = 
10 6 fibers and averaging over 10 3 samples both for uniform 
and Weibull distributions for the fiber breaking thresholds. 

Figure 2 displays the burst statistics for a fiber bun- 
dle with uniform and Weibull failure thresholds in Figs. 
2(a,c) and Figs. 2(6, d), respectively. In both cases, be- 
low the critical fraction of unbreakable fibers, a Cl the burst 
distributions D(A) do not change significantly, and even 
the cutoffs associated to the lack of numerical statistics do 
not change with a. For a uniform distribution of break- 
able fibers, the constitutive equation for a < a c (a) reduces 
to the one corresponding to a material composed only by 
weak fibers, with threshold values between zero and the 
upper bound <r 1 t f l ax = 1/ (1 — a). Hence, it follows that 
the entire failure process of the bundle, obtained at dif- 
ferent a values, remains the same until there are enough 
weak fibers in the system N w > N/2. Furthermore, for 
the parameter regime a < a c the avalanche statistics does 
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Fig. 3: The average size of the largest burst A ma x as a function 
of a for uniform (a, c) and Weibull (6, d) distributions. The 
vertical straight lines in the figures indicate the corresponding 
critical point a c - (c) and (d) show that approaching a c from 
above, A mal has a power law divergence as a function of a — a c . 
The value of the exponent is v = 1.56 ± 0.07 for both cases. 



not change. We obtain the typical power law distribution 
D(A) ~ A~ T with the exponent r = 5/2 (Figs. 2(a,c)). 
On the other hand, the parabolic shape of the constitutive 
behavior a (cr) also prevails for a > a c - In this regime, 
the system behaves as if the loading process was stopped 
before reaching the maximum of c (cr), due to the insuf- 
ficient number of breakable fibers N w < N/2 (compare 
to Fig. 1(a)). Consequently, the cutoff of the distribution 
D(A) in Fig. 2(a) decreases with increasing a. However, 
the exponent r keeps the same value as below a c , in agree- 
ment with our predictions and with Ref. [4]. 

We use the average size of the largest burst A max as the 
characteristic burst size of the system. It can be seen in 
Fig. 3(a) that for the uniform distribution below a c , the 
value of A max is constant, while it decreases rapidly when 
a surpasses a c . Figure 3(c) demonstrates that approach- 
ing a c from above the characteristic burst size shows a 
power law divergence A max ~ (a — a c )~ u . The value of 
the exponent v = 1.56 ± 0.07 was obtained numerically. 

In Figures 2(6,d), the burst statistics of a fiber bundle 
with a Weibull distribution of breaking thresholds is illus- 
trated. The Weibull parameters were set to m = 2 and 
A = 1, corresponding to a critical point a c ~ 0.3085. In 
this case, below the critical point a < a c the burst size 
distribution has a power law behavior D(A) ~ A~ r , with 
the exponent r = 5/2 (Fig. 2(6)). In this regime, a slight 
increase of the cutoff burst size appears when increasing a, 
however, the power law part of the distribution does not 
change. When a surpasses a c , the exponent of the power 
law regime of -D(A) suddenly switches to the lower value 
r = 9/4. The latter is in excellent agreement with our 



analytic predictions (see Fig. 2(6)). Moreover, we again 
find that the characteristic burst size A max diverges as a 
power law as we approach a c from above. The value of 
the critical exponent v is the same as for the uniform case 
(see Fig. 3(d)). We emphasize that the value of the expo- 
nent of the power law regime of D(A) remains constant 
r = 9/4 when changing a above a c . 

Using A rnax as a scaling variable, we introduce the scal- 
ing ansatz 

D(A) = A-^A/ALj (16) 

for the burst size distributions above the critical point a > 
a c . Here (3 and £ are scaling exponents, which have the 
relation j3 = r£ with r = 5/2 and r = 9/4 for the uniform 
and Weibull distributions, respectively. Figures 2(c) and 
(d) present the rescaled burst size distributions plotting 
f(A)A mra as a function of A/A^ aa: . The high quality 
data collapse is obtained with the parameters (5 — 3.25, 
£ = 1.25 and (3 = 2.52, £ = 1.12, for the uniform and 
Weibull distributions, which are consistent with the two 
different values of the r exponent. 

Discussion. — Our numerical and analytical calcula- 
tions revealed that the presence of unbreakable elements 
gives rise to a substantial change of the fracture process of 
disordered materials both on the micro- and macro-scales. 
Astonishingly we found a critical fraction of the break- 
able and unbreakable components where the exponent of 
the burst size distribution switches from the well known 
mean field exponent of FBM t = 5/2 to a significantly 
lower value r = 9/4. The transition is conditioned to 
disorder distributions where the macroscopic constitutive 
response of the system has a single maximum and an in- 
flexion point, implying a novel universality class of FBM. 
Despite we have considered only unbreakable fibers, our 
results will hold for a finite gap of two threshold distri- 
butions down to a certain critical value of the gap size. 
Below this critical value our model recovers the former 
work of Ref. [18]. Besides its theoretical importance, the 
problem has several implications for experimental stud- 
ies. New materials of high mechanical performance are 
often fabricated by mixing components with widely differ- 
ent properties. For instance, fiber reinforced composites 
are composed of strong fibers which are embedded in a 
carrier matrix. In this case at the breaking of weak ele- 
ments, the strong ones act as the unbreakable component 
of our model. 

Our detailed analytical and numerical study is restricted 
to the quasi-static limit of FBMs where the external load 
is incremented in a continuous manner. However, in lab- 
oratory experiments only finite discrete load steps can be 
realized. Recently, it has been demonstrated that for large 
enough load increments the statistics of bursts changes, i.e. 
under GLS conditions the exponent of the power law dis- 
tribution of burst sizes takes a higher value r = 3.0 [19,20]. 
When the system is a mixture of weak and strong fibers 
the effect of finite load steps depends on the type of dis- 
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order. For uniformly distributed failure thresholds it is 
straightforward to show that the burst exponent changes 
to t = 3.0 both below and above the critical fraction a c . 
Nevertheless, disorder distributions for which the constitu- 
tive curve has an inflexion point (e.g. Wcibull distribution) 
deserve a detailed study which will be presented elsewhere. 

A very interesting application of FBMs is to study the 
time dependent deformation and rupture of disordered 
materials under a constant external load (creep rupture). 
To understand damage enhanced creep processes the re- 
laxation dynamics of FBMs has recently been investigated 
in details [21-23]. It has been found that below the critical 
load <t c the system suffers only partial failure and relaxes 
to a stable state, while above a c macroscopic breaking oc- 
curs in a finite time. When approaching er c from cither 
side, the characteristic time scale (relaxation time and 
lifetime of the system) has a power law divergence with 
a universal exponent 1/2 [22,23]. Since the relaxation dy- 
namics of the system is determined by the functional form 
of the constitutive curve in the vicinity of the critical load, 
novel behaviour can be expected in the presence of strong 
fibers. The relaxation dynamics and creep rupture open 
up interesting possibilities for future applications of our 
model. 
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